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ABSTRACT 

We present the results of Chandra observations of six radio-loud quasars (RLQs) and one optically bright radio- 
quiet quasar (RQQ) at z ~ 4.1-4.4. These observations cover a representative sample of RLQs with moderate 
radio-loudness (R « 40-400), filling the X-ray observational gap between optically selected radio-quiet quasars 
(predominantly R <; 2-10) and the five known blazars at z > 4 {R « 800-27000), where R=f$ ghz//44oo A ( rest 
frame). We study the relationship between X-ray luminosity and radio-loudness for quasars at high redshift and 
constrain RLQ X-ray continuum emission and absorption. From a joint spectral fit of nine moderate-/? RLQs 
observed by Chandra, we find tentative evidence for absorption above the Galactic Nh, with a best-fit neutral in- 



trinsic column density of Nh = 2.4^'g x 10 /z cm~ z , consistent with earlier claims of increased absorption toward 
high-redshift RLQs. We also search for evidence of an enhanced jet-linked component in the X-ray emission due 
to the increased energy density of the cosmic microwave background (CMB) at high redshift, but we find neither 
spatial detections of X-ray jets nor a significant enhancement in the X-ray emission relative to comparable RLQs 
at low-to-moderate redshifts. Overall, the z ~ 4-5 RLQs have basic X-ray properties consistent with comparable 
RLQs in the local universe, suggesting that the accretion/jet mechanisms of these objects are similar as well. 



Subject headings: galaxies: active — galaxies: high-redshift - 
general 



- galaxies: jets — galaxies: nuclei — quasars: 



1. INTRODUCTION 

The study of the z > 4 universe has become a major theme 
in astronomy over the past decade. In particular, X-ray obser- 
vations of quasars at z > 4 reveal the physical conditions in the 
vicinities of the first massive black holes to form in the Uni- 
verse. Prior to 2000, there were only six published X-ray de- 
tections of quasars at z > 4 (see Kaspi, Brandt, & Schneider 
2000; hereafter KBS00). In the last four years, through the un- 
precedented sensitivities of Chandra and XMM-Newton as well 
as optical wide-area surveys such as the Sloan Digital Sky Sur- 
vey (SDSS; York et al. 2000), the number has risen to nearly 
one hundred. 

Equally important are wide-area radio surveys, such as Faint 
Images of the Radio Sky at Twenty Centimeters (FIRST; e.g., 
Becker et al. 1995) and the NRAO VLA Sky Survey (NVSS; 
Condon et al. 1998), to constrain the properties of radio cores, 
jets, and lobes. Radio-selected samples of radio-loud quasars 
(RLQs), although sampling a minority of the total quasar pop- 
ulation, are less prone to selection effects than optical sur- 
veys, as radio emission is not affected by absorption due to 
dust. The radio-loudness parameter, as defined by Kellermann 
et al. (1989), is given by R=fs ghz/ /4400 A ( res t frame), where 
quasars with R > 10 are RLQs and those with R < 10 are radio- 
quiet quasars (RQQs). 

Some of the first z > 4 quasars studied in X-rays were a small 
group of highly radio-loud blazars 4 (R « 800-27000) in which 
the X-ray radiation is probably dominated by a jet-linked com- 



ponent. These blazars represent a minority of the RLQ popu- 
lation, however, and so are not suitable for representative sta- 
tistical studies. Recently, Chandra's Advanced CCD Imaging 
Spectrometer (ACIS) has obtained X-ray detections of many 
RQQs (predominantly with R <J 2-10), which represent the ma- 
jority population of quasars in the local universe and so too, 
apparently, at high redshift (e.g., Schmidt et al. 1995; Stern et 
al. 2000). Through these studies (e.g., Brandt et al. 2002; Vig- 
nali et al. 2001, 2003a,b), the basic X-ray properties of z ~ 4-5 
RQQs have become reasonably well understood. Perhaps sur- 
prisingly, given the strong evolution of the quasar population 
(e.g., Boyle et al. 2000; Fan et al. 2001b), high-redshift RQQs 
have remarkably similar X-ray and broad-band properties to 
RQQs at low-to-moderate redshift (after controlling for lumi- 
nosity effects; e.g., Vignali, Brandt, & Schneider 2003, here- 
after VBS03), suggesting that the accretion mechanism does 
not dramatically evolve with time. Prior to the current paper, 
however, there were only four published X-ray detections of 
z > 4 RLQs with R « 10-1000 corresponding to typical RLQs 
in the local universe (excluding the blazar PMN 0525—3343 
which has R w 800). 

Here we report Chandra detections of six z ~ 4.1-4 .4 radio- 
selected RLQs and one RQQ (FIRST 0747+2739) which is ra- 
dio detected and optically bright. These objects have been se- 
lected as prime targets for X-ray observations based on their 
bright optical and radio fluxes, and they have been chosen to 
represent the majority population of RLQs. Four are FIRST- 
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selected objects published in Benn et al. (2002), two are from 
the Cosmic Lens All Sky Survey (CLASS; Myers et al. 2003) 
published in Snellen et al. (2001), and one is from the SDSS, 
first presented in Anderson et al. (2001). Five of the RLQs have 
published detections at both 1 .4 GHz and 5 GHz and have fiat 
radio spectra with a r ~ —0.4 to +0.3 (/ v °< V a ). In the case 
of FIRST 1309+5733, the upper limit on the 5 GHz flux den- 
sity gives a r < 0.06, so it is possible that this object is also a 
flat-spectrum RLQ. None of the quasars shows any significant 
structure in the FIRST radio images (these are the best radio im- 
ages available, with an angular resolution of m 5"). While the 
RQQ FIRST 0747+2739 is relatively radio-weak compared to 
our other objects, it is one of the optically brightest quasars at 
z > 4, and an overabundance of C IV absorption in the optical 
spectrum suggests that it may be a transition Broad Absorption 
Line (BAL) quasar (Richards et al. 2002). Given that RLQs 
are often substantially X-ray brighter than RQQs with compa- 
rable optical luminosity (e.g., Zamorani et al. 1981; Worrall et 
al. 1987), these objects are prime targets for Chandra snapshot 
(e.g., ~ 4-10 ks) observations. Indeed, all seven of our sources 
are detected with a minimum of « 20 counts. 
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FIG. 1. — Radio loudness versus redshift for X-ray detected quasars at 
Z > 4. The large symbols represent the moderate-S RLQs (triangles), blazars 
(stars), and RQQs (circles). The filled symbols represent the quasars in the 
current sample. Arrows show upper limits for radio loudness at the ks 3a level 
based on the 1.4 GHz radio flux density. Note that, excluding upper limits, 
the present sample more than doubles the number of objects with moderate 
radio-loudness (R 10-1000), filling the X-ray observational gap between the 
RQQs and the blazars. 



representative sample for reliable statistical evaluations of z > 4 
RLQ X-ray properties. 

Recent work (e.g., Cirasuolo et al. 2003; Ivezic et al. 2002, 
2004) has sparked new controversy over the existence of bi- 
modality in the radio-loudness distribution of quasars, sug- 
gested first by Strittmatter et al. (1980). Current results suggest 
that, although quasars cover the full range of radio-loudness, 
the underlying distribution does consist of two regimes with a 
local minimum in between. As reported by Ivezic et al. (2002) 
for a sample of 3225 SDSS quasars matched to FIRST sources, 
8% ± 1% of all quasars with ;'* < 18.5 are radio-loud, appar- 
ently independent of optical luminosity or redshift, in the range 
z < 2. Figure 2 shows the radio-loudness distribution for z > 4 
quasars with radio detections or upper limits in FIRST, with the 
best-fit distribution from Ivezic et al. (2004) overlaid (dashed 
curve). While the blazars lie on the extreme radio-loud tail of 
the distribution, the RLQs from this work fall near the peak of 
the radio-loud regime, likely representing the majority popula- 
tion of z > 4 RLQs. By integrating the best-fit curve directly, 
we find that our sample of RLQs with R = 40-400 is represen- 
tative of A 50% of the total RLQ population, while the blazars 
with R > 1000 represent < 5% of all RLQs. 
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This brings the total sample of z > 4 X-ray detected RLQs 
to 15, more than doubling the number of moderately radio- 
loud quasars and filling the gap (as shown in Figure 1) between 
the many X-ray observations of z > 4 RQQs and the extreme 
blazars. 5 To date there are « 50 known z > 4 RLQs; thus the 
objects in our sample comprise w 10% of those known. Our tar- 
gets have been drawn from surveys covering a significant frac- 
tion of the sky (6400 deg 2 , Snellen et al. 2001 ; 4030 deg 2 , Benn 
et al. 2002) and so are likely to be representative of z > 4 RLQs 
as a whole, constituting the first set large enough to serve as a 

5 A regularly updated list of z, > 4 quasars with X-ray detections can be found at http 
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FIG. 2. — Radio-loudness distribution for known z > 4 quasars with de- 
tections or upper limits in the FIRST survey. The solid histogram shows all 
z > 4 quasars with radio detections, and the dotted histogram shows all those 
with radio upper limits. The dashed curve is adapted from the radio-loudness 
distribution in Ivezic et al. (2004), determined from a sample of ~ 10,000 ob- 
jects detected by FIRST and SDSS (with arbitrary normalization). The objects 
in the current sample (hatched shading) represent the majority population of 
RLQs, while the blazars (solid shading) lie on the extreme radio-loud tail. 

Throughout this paper we adopt Hq=10 km s _1 Mpc -1 in 
a A-cosmology with il M =0.3 and i2 A =0.7 (e.g., Spergel et al. 
2003). 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Chandra Observations and Basic Data Reduction 

The X-ray observations were acquired by Chandra dur- 
ing Cycle 4, using the Advanced CCD Imaging Spectrometer 
(ACIS; Garmire et al. 2003) with the S3 CCD at the aimpoint. 

//www.astro.psu.edu/users/niel/papers/highz-xray-detected.dat 
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Table 1 

Chandra OBSERVATION LOG 



Object 




FIRST 


FIRST 


ARad-X a 


X-ray 


Exp. Time b 




Name 


z 




«2000 


82000 


(arcsec) 


Obs. Date 


(ks) 


Ref. 


FIRST 0725+3705 


4.33 


07 


25 18.26 


+37 05 18.3 


0.2 


2002 Dec 17 


4.89 


(1) 


FIRST 0747+2739 


4.11 


07 


47 11.20 


+27 39 04.0 


0.9 


2002 Dec 03 


4.96 


(1,3) 


SDSS 0839+5112 


4.39 


08 


39 46.20 


+51 12 02.9 


0.3 


2003 Jan 23 


4.90 


(2,4) 


CLASS J0918+0637 


4.19 


09 


18 24.39 


+06 36 53.3 


0.2 


2002 Dec 10 


4.90 


(2) 


FIRST 1309+5733 


4.27 


13 


09 40.69 


+57 33 10.0 


0.4 


2003 Jul 07 


4.64 


(1) 


CLASS J1325+1 123 


4.40 


13 


25 12.48 


+ 11 23 30.0 


0.5 


2003 Mar 02 


4.70 


(2) 


FIRST 1423+2241 


4.32 


14 


23 08.24 


+22 41 58.2 


0.2 


2002 Dec 02 


4.70 


(1) 



NOTE. — Radio and optical identifications of the quasars presented here can be found in the papers cited in the reference 
column. Radio positions are from the FIRST catalog (Becker et al. 1995). Units of right ascension are hours, minutes, 
and seconds, and units of declination are degrees, arcminutes, and arcseconds. 

a Distance between the radio and X-ray positions. b The Chandra exposure time has been corrected for detector dead 
time. 



REFERENCES. — (1) Benn et al. 2002; (2) Snellen et al. 2001; (3) Richards et al. 2002; (4) Anderson et al. 2001. 



The observation log is reported in Table 1 . Faint mode was used 
for the event telemetry format, and the ASCA grade 0, 2, 3, 4, 
and 6 events were used in the analysis. With one exception, the 
observations are free from background flares. A flare occurred 
near the end of the observation of SDSS 0839+51 12, resulting 
in an increased background count rate (by a factor of « 10) for 
w 10 minutes. Due to the superb imaging capabilities of Chan- 
dra, this flare does not have a substantial effect on our data; the 
count rates used in deriving physical quantities have been back- 
ground subtracted, and, especially since SDSS 0839+5112 is 
the X-ray brightest source of our sample, the additional back- 
ground does not significantly affect our spectral analysis. 

Source detection was carried out with WAVDETECT (Free- 
man et al. 2002). For each image, we calculated wavelet trans- 
forms (using a Mexican hat kernel) with wavelet scale sizes 1, 
1.4, 2, 2.8, and 4 pixels. Those peaks whose probability of be- 
ing false were less than the threshold of 10~ 5 were declared 
real. This threshold is entirely appropriate when, as is the case 
here, the source position of interest is specified a priori. The 
Chandra positions of the detected quasars lie within 0.2-0.9"of 
their FIRST radio positions (see Table 1), consistent with the 
expected positional error. Source detection was performed in 
four energy ranges: the ultrasoft band (0.3-0.5 keV) , the soft 
band (0.5-2 keV), the hard band (2-8 keV) and the full band 
(0.5-8 keV). In the redshift range z«4. 1^4.4 for these objects, 
the full band corresponds to the « 2.5^-3 keV rest-frame band. 

The WAVDETECT X-ray photometry is shown in Table 2. We 
have checked these results against manual aperture photome- 
try and found good agreement between the two methods. Ta- 
ble 2 also lists the hardness ratios [defined as (H — S)/(H + S), 
where S is the soft-band (0.5-2 keV) counts and H is the hard- 
band (2-8 keV) counts], the band ratios (H/S), and the effec- 
tive X-ray power-law photon indices (r) derived from these 
band ratios assuming Galactic absorption. This derivation has 
been performed using the Chandra X-ray Center Portable, In- 
teractive, Multi-Mission Simulator (PIMMS; Mukai 2002) soft- 
ware, and we have applied a time-dependent correction for the 
quantum-efficiency degradation of Chandra ACIS at low ener- 
gies (see §3 for details). All of the quasars have effective pho- 
ton indices consistent, within the significant errors, with previ- 
ous studies of RLQs at z ~ 0-2 (r w 1.4-1.9; e.g., Reeves & 
Turner 2000). We also performed basic spectral analyses for all 
of the sources using the Cash statistic (Cash 1979) with XSPEC 



(Version 1 1.2.0; Arnaud et al. 1996) in the 0.5-8 keV band. We 
found general agreement between the two methods. Further 
spectral analysis is presented in §3.2. 

Figure 3 shows the soft-band images (chosen to give the best 
signal-to-noise ratio) of the six RLQs, both as raw 0.492" pix- 
els and adaptively smoothed at the 3o level using the algorithm 
of Ebeling, White, & Rangarajan (2004). 

2.2. Companion X-ray Sources 

Powerful X-ray jets have been observed emanating from 
many RLQs, 6 and their presence has been inferred in high- 
redshift RLQs (particularly the blazars, in which the jet is pre- 
sumably oriented close to the line of sight) from high lumi- 
nosity and variability in both the X-ray and radio bands (e.g., 
Fabian et al. 1997, 1999). At high redshift such investigations 
are particularly interesting, since if a significant component of 
the X-ray radiation is a consequence of inverse Compton (IC) 
scattering of cosmic microwave background (CMB) photons, 
then the X-ray emission may be significantly enhanced (see 
§3.3 for details). Additionally, in the favored cosmological 
models, the angular size increases with redshift beyond z ~ 1 .6 
such that the sub-arcsecond spatial resolution of Chandra is 
sufficient to resolve kpc-scale jets at all redshifts. Recently, 
X-ray jets have been discovered associated with the z = 3.8 ra- 
dio galaxy 4C41.17 (Scharf et al. 2003) and with the z = 4.3 
blazar GB 1508+5714 (Siemiginowska et al. 2003; Yuan et al. 
2003). Schwartz (2004) also presents convincing evidence for 
the detection of an IC/CMB enhanced jet associated with the 
z = 4.01 blazar GB 1713+2148 in a ROSAT HRI observation. 

We have searched for possible companions and jets over a 
region of « 100" x 100" centered on each quasar (the same re- 
gions shown in the smoothed images in Fig. 3). Atz ~ 4.1-4.4, 
100" corresponds to a linear scale of w 700 kpc. In the seven 
fields, we found a total of three objects in the soft band, which 
can be seen in Fig. 3. None of these is sufficiently X-ray bright 
and close to our targets to contaminate their X-ray emission. 
The three serendipitous detections visible in Fig. 3 all have 
B-band counterparts in the second Palomar Optical Sky Sur- 
vey (POSS II) plates, 7 which probably precludes them from 
being companion objects to the quasars at high redshift. In ad- 
dition, we detected a faint (3 counts) X-ray source in the hard 
(2-8 keV) band « 34" to the southeast of SDSS 0839+5112, 
which appears to have faint counterparts in both the POSS II 



6 A list of X-ray detected quasar jets is available at http://hea-www.harvard.edu/XJET/ 

7 Available at http://www.nofs.navy.mil/data/fchpix/ 
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Table 2 

X-ray Counts, Hardness Ratios, Band Ratios, and Effective Photon Indices 

X-ray Counts" 



Object [0.3-0.5 keV] [0.5-2 keV] [2-8 keV] [0.5-8 keV] Hardness Ratio b Band Ratio b T c 



FIRST 0725+3705 
FIRST 0747+2739 
FIRST 1309+5733 
FIRST 1423+2241 
CLASS 10918+0636 
CLASS 11325+1123 
SDSS 0839+5112 


<3.0 

<4.8 

<4.7 

<3.0 

5.0™ 

<4.5 

2.0+|j 


21-2^ 
16.8j*g 
24 8+ 61 
16.8+2? 
41.6+12 
23.7121 
55.0+H 


^■ y -2.1 

5 0+ 34 

6 9+ 3 ' 8 
ay -2.6 

2 Q+ 2 - 9 
13.7+^ 

20.612s 


26.1+gf 
21 1 +5J 
31.7+j; 7 6 
19.7+22 

55.4l?2 
30.6+J5 
75 4+ 9:7 


-0 62+ uw 
uoz -0.14 

-0 54+ 22 

_ 057 +0.18 
u --"-0.14 

-0.79+°^ 

-o. 5 o+° : ; 3 

u -+ J -o.io 


0.23+11? 

0.30+°;" 

0.28+°;}« 

12+°' 16 
u - lz -0.08 

33+ 012 
u - J -'-o.io 

0.29+JJ;! 7 
0.38+» : li 


9 ,+0.fc 
za -0.5 
1.8+»j 

2 6+ 10 
z -°-0.8 

1.8+0.3 

1 R+°- 5 
1 - S -0.4 

1 6+ 03 


a Errors on the X-ray counts were computed according 
calculated using Poisson statistics. The upper limits are 


to Tables 1 and 2 of Gehrels (1986) and correspond to the 10 level; these were 
at the 95% confidence level and were computed according to Kraft, Burrows, & 



Nousek(1991). 

b Errors on the hardness ratios [defined as (H — S)/(H + S), where S is the soft-band (0.5-2 keV) counts and H is the hard-band 
(2-8 keV) counts], the band ratios (H/S), and the effective photon indices are at the ks la level and have been computed following the 
"numerical method" described in § 1.7.3 of Lyons (1991). This avoids the failure of the standard approximate variance formula when 
the number of counts is small (see § 2.4.5 of Eadie et al. 1971). 

c The effective photon indices (and errors) have been computed from the band ratios and their respective errors, using the PIMMS 
software. We have applied a time-dependent correction which accounts for the quantum-efficiency degradation of Chandra ACIS at low 
energies (see §3). 



B-band and FIRST images (no B-band counterpart). However, 
since no counts were detected at its position in the soft band and 
jets typically have fairly soft X-ray spectra (r w 1.6-2.3; Har- 
ris & Krawczynski 2002), it appears unlikely that the source is 
a jet from SDSS 0839+5112. 

As a further observational check, we compared the 
0.5-2 keV cumulative number counts of the sources detected 
in the quasar fields (the targeted objects were, of course, ex- 
cluded) with those presented by Moretti et al. (2003) using a 
compilation of six different surveys. The surface density of 
X-ray sources at our observed flux limit (w2x 10~ 15 erg cm~ 2 
s~ ! ) is N(> S) = 590^320 deg~ 2 , consistent within the errors 
with the results (490 deg~ 2 at this flux limit) of Moretti et al. 
(2003). Thus we find no evidence for an excess of companion 
objects associated with our quasars. 

2.3. X-ray Extension of the Quasars 

To constrain the presence of either gravitational lensing (e.g., 
Wyithe & Loeb 2002; Comerford, Haiman, & Schaye 2002) or 
X-ray jets close to the sources, we performed an analysis of 
X-ray extension for all the quasars in the present sample. For 
each quasar, we created a point spread function (PSF) at « 1.5 
keV at the source position (using the CIAO tool MKPSF) nor- 
malized by the observed number of source counts. Then we 
compared the radial profile of this PSF with that of the source, 
finding agreement for each object. To identify a putative jet 
by this method, we require a minimum of m 3 nearby counts 
offset by « 2" from the core, corresponding to a distance of 
« 15 kpc at z = 4.1-4.4. We can exclude significant X-ray ex- 
tension above this scale for the objects in the present sample. 

2.4. X-ray Variability 

Given tentative reports of increasing quasar X-ray variabil- 
ity with redshift (in the sense that quasars of the same X-ray 
luminosity are more variable at z > 2; Manners, Almaini, & 
Lawrence 2002), we have searched for X-ray variability by an- 
alyzing the photon arrival times in the 0.5-8 keV band using 
the Kolmogorov-Smirnov (KS) test. While no significant vari- 
ability was detected (the KS probabilities indicate that a con- 
stant hypothesis can only be rejected at the w 22-92% level), 

8 This list is available from http://www.astro.psu.edu/users/niel/papers/papers.html 



we note that the analysis of X-ray variability for the quasars 
studied here is limited by both the short X-ray exposures in the 
rest frame (« 15 min; see Table 1) and the low X-ray fluxes of 
most of our quasars. 

3. X-RAY ANALYSIS 

The principal X-ray, optical, and radio properties of the fif- 
teen z > 4 RLQs with X-ray detections (including those de- 
tected prior to this work) are given in Table 3. A description is 
as follows: 

Column (1 ). — The name of the source. 

Column (2). — The redshift of the source. 

Column (3). — The Galactic column density (from Dickey & 

Lockman 1 990) in units of 1 20 cm" 2 . 

Column (4). — The monochromatic rest-frame ABi45on +z ) 
magnitude (defined in §3b of Schneider et al. 1989), cor- 
rected for Galactic extinction. When spectroscopically deter- 
mined values are not available, AB 1450 ( 1+Z ) magnitudes have 
been derived from the SDSS /-magnitudes assuming the em- 
pirical relationship AB i +z \ = i — 0.2. For objects with- 
out SDSS detections, we used B-band magnitudes from the 
APM Catalog (McMahon et al. 2002) and the relationship 
AB 1450(1+Z ) =R- (0.648)z + 3.10. We have compiled a list of 
all z > 4 quasars with similar photometry; this is used in the 
analysis throughout this paper. 8 All of the above relationships 
provide reliable AZ?i45o(i +z ) estimates (within w 0.1-0.2 mags) 
in the redshift range under consideration. 
Columns (5) and (6). — The 2500 A rest-frame flux density 
and luminosity. These were computed from the AB 1450 ( 1+Z ) 
magnitude assuming an optical power-law slope of a = —0.5 
(/ v cx v a ; Vanden Berk et al. 2001). The 2500 A rest-frame 
flux densities and luminosities are increased by « 15% for an 
optical power-law slope of a = —0.79 (e.g., Fan et al. 2001a) 
as in Vignali et al. (2001, 2003a,b). 

Column (7). — The absolute B-band magnitude computed as- 
suming a = —0.5. If a = —0.79 is adopted for the extrapo- 
lation, the absolute B-band magnitudes are brighter by w 0.35 
mag. 

Columns (8) and (9). — The observed count rate in the 
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FIG. 3. — Chandra 0.5-2 keV 2.5-10 keV rest frame) images of the six RLQs presented in this paper. In each panel, the horizontal axis shows the Right 
Ascension, and the vertical axis shows the Declination (both in J2000 coordinates). North is up, and East to the left. For each pair, the upper images are 97" x 97" 
and have been adaptively smoothed at the 3c level [except for FIRST 1309+5733 (2c level)] using the algorithm of Ebeling et al. (2004). The lower images are 
10" x 10", showing raw 0.492" pixels. Crosses mark the optical positions of the quasars. 



0.5-2 keV band (for Chandra observations) and the corre- 
sponding flux (fx), corrected for Galactic absorption. This flux 
has been calculated using PIMMS and a power-law model with 
r = 1.6, as derived for samples of z ~ 0-2 RLQs (e.g., Reeves 
& Turner 2000); see also §3.2. Changes of the photon index 



in the range r= 1.4-1.9 lead to only a few percent change in 
the measured X-ray flux. The soft X-ray flux derived from the 
0.5-2 keV counts is generally similar (to within w 5-15%) to 
that derived using the full-band (0.5-8 keV) counts. 

The X-ray fluxes reported and used in this paper have been 
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FIG. 4. — X-ray - optical correlations for z > 4 RQQs (circles), RLQs (triangles), blazars (stars) and X-ray upper limits (downward-pointing arrows, at the 
ss 3a confidence level). The new Chandra observations presented in this paper are marked with filled symbols, (a) Observed-frame, Galactic absorption-corrected 
0.5-2 keV flux versus AB 1450 j 1+z ) magnitude for z > 4 quasars. The dashed lines show the a ox = -1.5 and a ox = —1-8 loci at z = 4.3 (the average redshift of the 
present sample). The dotted line shows the best-fit correlation reported in §4 [Equation (2)] of Vignali et al. (2003b) for z > 4 optically selected RQQs. (b) X-ray 
(2 keV) versus UV (2500 A) luminosity density for z, > 4 quasars. The dotted line indicates the best-fit relationship for 137 SDSS RQQs in the 0.16-6.28 redshift 
range from VBS03. Clearly, the RLQs lie significantly above the best-fit line, typically being X-ray brighter by a factor of ss 2 than RQQs at comparable UV 
luminosities. 



corrected for the ACIS quantum-efficiency degradation at low 
energy using the method described in §3.2. The 0.5-2 keV flux 
correction is w 25% for the Chandra observations presented 
here. 

Columns (10) and (11). — The rest-frame 2 keV flux density 
and luminosity, computed assuming F = 1.6. 
Column (12). — The 2-10 keV rest-frame luminosity, cor- 
rected for Galactic absorption. 

Column (13). — The optical-to-X-ray power-law slope, oc ox , 
defined as 

„ _ log(/2kev//2500 A) n n 

log(V2kev/V 2 500 A) 

where fa k eV and /2500 A WQ the rest-frame flux densities at 
2 keV and 2500 A, respectively. The « la errors on a ox have 
been computed following the "numerical method" described in 
§ 1.7.3 of Lyons (1991). Both the statistical uncertainties on 
the X-ray count rates and the effects of possible changes in the 
X-ray (T w 1.4-1.9) and optical (a « -0.5 to -0.9; Schneider 
et al. 2001) continuum shapes have been taken into account. 
Our choice of F = 1.6 and a = —0.5 accounts for the slight dif- 
ferences in fasoo, /2kev, and a ox values as compared to those in 
Vignali et al. (2001, 2003a,b). 

Column (14). — The radio power-law slope (fa v a ) between 
1 .4 GHz and 5 GHz (observed frame). The 1 .4 GHz flux density 
is from FIRST or NVSS (Condon et al. 1998), and the 5 GHz 
flux density is from the GB6 (8 > 0; Gregory et al. 1996) or 
PMN (8 < 0; Griffith & Wright 1993) surveys. 
Column (15). — The radio-loudness parameter defined by 
Kellermann et al. (1989) as R — fa ghz/ /4400 A ( rest frame). 
The rest-frame 5 GHz flux density was computed from the 
FIRST or NVSS observed 1.4 GHz flux density and the radio 
power-law slope given in column 14. Where a 5 GHz detection 



was not achieved, we have assumed a flat spectrum of a r = to 
be consistent with the other z > 4 RLQs. The rest-frame 4400 A 
flux density was computed from the AB 14 5 ( 1+Z ) magnitude as- 
suming an optical power-law slope of a = —0.5. 
Column (16). — Recent X-ray references. 

3.1. X-ray - Optical Correlations 

With nearly one hundred published X-ray detections of 
quasars at z > 4, the typical X-ray fluxes of z « 4-5 quasars 
are now well defined. Plots of X-ray flux versus AB 14 5 (i +z ) 
magnitude are often constructed for distant quasars and are 
useful for planning future observations, as they reflect directly 
observable quantities. Figure 4a is such a plot, depicting the 
observed-frame 0.5-2 keV flux versus AB 145 o( 1+z - ) for quasars 
at z > 4 with X-ray detections and upper limits. A strong cor- 
relation exists between these quantities, and the best-fit rela- 
tionship derived in Vignali et al. (2003b) for z > 4 RQQs is 
shown as a dotted line. It is immediately apparent that the 
RLQs in this sample (depicted as filled triangles) are generally 
brighter in X-rays than RQQs (circles) with comparable optical 
magnitudes, all lying above the best-fit relationship. The RQQ 
FIRST 0747+2739 (filled circle) lies below the best-fit. 

Extending much earlier work, VBS03 found that the X-ray 
and broad-band properties of RQQs do not significantly evolve 
with redshift. VBS03 found that, independent of redshift in 
the range z ~ 0.2-6.3, the X-ray luminosities (Lx) of RQQs are 
highly correlated (9.3a in VBS03) with their ultraviolet (UV) 
luminosities (Ltjv)- These quantities are plotted for the X-ray 
detected z > 4 quasars in Figure 4b, along with the best-fit re- 
lationship for RQQs from VBS03 (adapted for an optical slope 
of a = —0.5). Clearly, the RLQs (triangles) lie significantly 
above the best-fit line for RQQs (circles), with systematically 
higher X-ray luminosities by a factor of « 2. These are the 




FIG. 5 . — Combined spectrum (in the observed frame and used only for presentation purposes) for the nine Chandra observations of z > 4 RLQs with moderate-R, 
fitted with a power-law model and Galactic absorption (see §3.2 for details). The 68, 90, and 99% confidence regions for the photon index and intrinsic column 
density obtained using the Cash statistic are shown in the inset. 



same trends exhibited by samples of RLQs at low-to-moderate 
redshift, in which the additional X-ray emission is attributed 
to a jet-linked component. In general, the excess X-ray emis- 
sion tends to correlate with radio-loudness (see, e.g., Figure 6), 
lending additional support to this interpretation. We further in- 
vestigate the implications of these trends and the nature of the 
jet-linked X-ray emission in §3.3 

3.2. Joint Spectral Fitting 

To provide tighter constraints on the average X-ray spectral 
properties of z > 4 RLQs, we have performed a joint spectral 
analysis of the nine moderate-/? (R w 4CM-00) RLQs observed 
by Chandra; these have an average redshift of (z) =4.43. 
Three of these objects have been previously published: 
PSS 0121+0347, Vignali et al. (2003a); SDSS 0210-0018, 
Vignali et al. (2001); and SDSS 0913+5919, Vignali et al. 
(2003b). Source counts have been extracted from 2" radius cir- 
cular apertures centered on the X-ray position of each quasar. 
The background has been taken from annuli centered on the 
sources, avoiding the presence of nearby faint X-ray sources. 
Only one object is characterized by low counting statistics 
(SDSS 0913+5919, which has 6 counts); all of the others have 
w 20-80 counts in the 0.5-8 keV band. The sample does not 

9 See http://www.astro.psu.edu/users/chartas/xcontdir/xcont.html 



appear to be biased by the presence of a few objects with atyp- 
ically high signal-to-noise ratios (see Table 2), and the removal 
of any one of these objects does not produce significantly dif- 
ferent results. 

We have corrected each spectrum for the quantum-efficiency 
degradation of ACIS at low energies using the correction sup- 
plied via the Chandra Calibration Database (CALDB), Version 
2.26. While the degradation is not severe above w 0.7 keV, it 
could significantly affect estimates of X-ray absorption if not 
properly corrected. We have also repeated our analysis using 
the ACISABS correction, 9 based on data from the external cal- 
ibration source aboard Chandra (the CALDB correction is de- 
rived from gratings data). Due to differences in the extrapo- 
lation of the time-dependent quantum efficiency between the 
models (see Marshall et al. 2004), the ACISABS correction leads 
to fits with a slightly lower absorption column and photon in- 
dex, but the results are generally consistent. 

The joint spectral fitting of the nine unbinned quasar spectra 
was carried out using the Cash statistic (Cash 1979), which is 
well suited to low-count sources (e.g., Nousek & Shue 1989). 
Since the Cash statistic is applied to the unbinned data, it al- 
lows us to retain all available spectral information. We can 
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also assign to each source its own Galactic absorption and red- 
shift. The resulting joint spectrum (in the « 2-40 keV rest- 
frame band), while derived from only « 350 counts, is reason- 
ably good due to the extremely low background of Chandra 
(there are only w 5 background counts in all source cells). Er- 
rors in the following analysis are quoted at 90% confidence for 
one interesting parameter. 

Fitting a power-law model with Galactic absorption using 
XSPEC, we find a good fit (shown in Figure 5) with an av- 
erage photon index of Y — 1.65 ±0.15. The fitted photon in- 
dex is consistent with that derived using the band ratios for 
the six RLQs in Table 2, from which we obtain a weighted 
mean of T= 1.77 ±0.16. It is also consistent with previous 
studies of z ~ 0—4 RLQs (e.g., Brinkmann et al. 1997; Reeves 
& Turner 2000) which have found that core-dominated, flat- 
spectrum RLQs typically have F w 1.4-1. 9. 

To search for evidence of intrinsic X-ray absorption in the 
above sample of nine moderate-/? RLQs, we added to the joint 
spectral model a redshifted neutral absorption component (the 
XSPEC model zphabs) with the redshift set to that of each 
source. With this model, we find tentative evidence ( £ 90% 
confidence) for absorption in excess of the Galactic Nh, with 



a best-fit value of N H = 2.4+ 2 '° x 10 ZZ cm~ z (errors at 90% 
confidence for one interesting parameter). This fit produces 
a somewhat higher photon index, T = 1.80^{j 25, which is still 
consistent with the typical range for RLQs from above. The 
confidence contours in the Nh-F plane are shown in Figure 5. 

From a large sample of ASCA observations, Reeves & Turner 
(2000) confirmed earlier suggestions by Cappi et al. (1997) 
and Fiore et al. (1998) of a significant Nh-z correlation for 
RLQs, finding that z ~ 2-4 RLQs have typical column den- 
sities of a few xl0 22 cm~ 2 (see their Figure 9). Also, for 
two of the z > 4 blazars, strong intrinsic X-ray absorption 
(10 22 -10 23 cm~ 2 ) has been suggested by spectral flattening at 
low energies: PMN 0525-3344 (Worsley et al. 2004a) and 
GB 1428+4217 (e.g., Fabian et al. 2001; Worsley et al. 2004b). 
Our findings are generally consistent with these claims and ex- 
tend this result to the majority of RLQs with moderate radio- 
loudness (the average radio-loudness of the objects in the above 
samples is (R) « 10000, while that for the present sample is 
(R) « 200). 

3.3. The Nature of the Jet-Linked Component 

The enhanced radio and X-ray emission of RLQs in com- 
parison to RQQs is usually attributed to a relativistic jet-linked 
component. For core-dominated, flat-spectrum RLQs, the ma- 
jority of the X-ray emission is presumably from the small-scale 
jet close to the AGN, which, aided by Doppler boosting, ac- 
counts for the increased X-ray luminosities and flatter X-ray 
spectra of RLQs as a whole (e.g., Worrall et al. 1987; Wilkes 
& Elvis 1987). For more extended (kpc scale) X-ray jets, there 
has been considerable debate over the process responsible for 
the X-ray emission. In several cases, low optical luminosi- 
ties, which preclude pure synchrotron emission from a single 
electron population, have led to synchrotron models with more 
complicated electron energy distributions (e.g., Stawarz et al. 
2004) or models based on some form of IC scattering. In the 
case of IC emission, individual observations of X-ray jets lead 
to differing conclusions as to the primary source of IC seed 
photons (see Harris & Krawczynski 2002 for a review). The 
obvious choices are the synchrotron photons themselves (syn- 
chrotron self-Compton - SSC), UV photons from the accretion 



disk, or cosmic microwave background photons (IC/CMB); the 
relative importance of these sources may depend on factors 
such as the quasar's environment, orientation, and the bulk rel- 
ativistic motion of the jet. Both SSC and IC from UV pho- 
tons produced in the AGN presumably dominate close to the 
core of the quasar, while IC/CMB may dominate at large dis- 
tances. At high redshift IC/CMB is particularly relevant since, 
as pointed out originally by Rees & Setti ( 1 968) and recently by 
Schwartz (2002), the cosmological decrease in surface bright- 
ness as (1 +z)~ 4 is compensated for by a (1 +z) 4 increase in 
the energy density of the CMB. Above z ~ 3-4, X-ray jets dom- 
inated by the IC/CMB process could even outshine their parent 
quasars (see Figure 2 of Schwartz 2002). 
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FIG. 6. — Broad-band spectral index a ox versus radio-loudness (R) for sam- 
ples of flat-spectrum RLQs at differing redshifts. The small circles represent 
the RLQs from the radio-selected sample of Brinkmann et al. (1997), at red- 
shifts of z ~ 0.2-2 (filled circles) and z ~ 2-4 (open circles). The large sym- 
bols represent z > 4 RLQs with X-ray detections and upper limits (downward 
pointing arrows). We see no significant evidence for X-ray enhancement of 
RLQs at high redshift, as both the z > 4 moderate-R RLQs (triangles) and the 
blazars (stars) fall well within the distribution of objects at lower redshift. The 
new Chandra observations from this paper are plotted as filled triangles. 



As established in §3.1, our z > 4 RLQs are w 2 times X-ray 
brighter than RQQs with comparable optical luminosity, so the 
X-rays we have detected appear to be largely from a jet-linked 
component. We do not, however, observe any evidence for spa- 
tially extended emission in the Chandra images (see §2.3), nor 
do the FIRST radio images show any obvious spatial structure 
(note that the sub-arcsecond angular resolution of Chandra is 
better than the « 5" resolution of the FIRST survey). In Fig- 
ure 6 we have plotted the broad-band spectral index a 

ox versus 

radio loudness for our observations of z > 4 RLQs (triangles) 
along with the Brinkmann et al. (1997) sample of radio-selected 
flat-spectrum RLQs which spans the z < 2 (solid circles) and 
z = 2-4 (open circles) redshift ranges. Although the number 
of moderately radio-loud RLQs in the Brinkmann et al. (1997) 
sample is limited, it is clear that the z > 4 RLQs are compa- 
rable in their X-ray and broad-band properties to those at low- 
to-moderate redshift. The z > 4 blazars (star symbols) also lie 
within the observed range of these parameters for comparable 
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low-redshift objects. We have performed a similar comparison 
of our z > 4 RLQs to the low-to-moderate redshift sample of 
flat-spectrum RLQs presented in Fig. 4 of Fabian et al. (1999) 
and reach the same conclusion. 

Figure 6 has implications for jet X-ray emission mechanisms 
in RLQs. In particular, the observations provide evidence that 
the IC/CMB process does not dominate the total X-ray pro- 
duction of these objects; if it did, then the the arguments of 
Schwartz (2002) predict that these jets at z > 4 should be as 
bright or brighter in X-rays than the quasar cores. Such an 
enhancement would lead to a systematic flattening of oc ox by 
« 0.1-0.2 relative to comparable objects at low redshift, which 
we do not observe. Rather, it seems that the X-rays are probably 
made primarily on small scales where the quasar photon field 
dominates the CMB photon field. This is consistent with ex- 
pectations for low-redshift flat-spectrum RLQs, where on small 
scales SSC and external IC scattering from AGN seed photons 
are thought to dominate. 

Still, the question remains: if we believe that much of the 
X-ray emission from these objects is associated with a jet- 
linked component on small scales, then why do we not see 
the extended IC/CMB jet, which should be significantly en- 
hanced by the increased energy density of the CMB at high red- 
shift? Schwartz (2002) uses the jets from the quasars 3C 273 
(z = 0.158) and PKS 0637-752 (z = 0.654) as examples of 
IC/CMB jets (at projected distances of roughly 25 kpc and 
75 kpc, respectively) that should outshine their parent quasars 
above z ~ 4. Both of these quasars have flat (core) radio spectra 
similar to those of our objects, and the predicted Chandra count 
rates of the cores as they would appear at z ~ 4 (see Figure 2 of 
Schwartz 2002) of a fewx 10~ 3 cts s _1 are consistent with our 
observations. However, 3C 273 and PKS 0637—752 are more 
radio loud than our quasars by nearly an order of magnitude, 
with R w 1500 and 2200, and so it may be that differences in 
either orientation or intrinsic jet power may account for our ob- 
servations. 10 But even in the case of the recently discovered 
X-ray jet from the z = 4.3 blazar GB 1508+5714, where the 
IC/CMB process has been invoked (Siemiginowska et al. 2003; 
Yuan et al. 2003; Cheung 2004), we do not see such a dramatic 
enhancement. This jet only accounts for w 3% of the quasar 
luminosity and was discovered in a Chandra observation with 
a lengthy (89 ks) exposure time. To test our sensitivity to de- 
tecting such jets, we randomly selected 800 s segments of the 
GB 1508+5714 data, chosen to have « 40 counts in the full 
0.5-8 keV band, matching the average number of counts in our 
w5 ks Chandra snapshot observations. We analyzed these im- 
ages following the procedure of §2.3 and did not find clear evi- 
dence for extent in any case. Inspecting the raw and adaptively 
smoothed images by eye, only in w 10-20% of the fields could 
we observe hints of the jet. This is not surprising, as at w 3% 
of the quasar luminosity we expect only 1-2 counts from the jet 
in each field. So if our moderate-/? RLQs have jets similar in 
nature to that of GB 1508+5714 (perhaps at a larger viewing 
angle to account for the difference in radio-loudness) then we 
would not expect to detect extended jets from these snapshot 
images. 

Overall, from this preliminary sample, it appears that the 
small-scale jet components of flat-spectrum RLQs do not 
change significantly between redshifts and 5. This result is 

10 The lower radio loudnesses of our objects may indicate that they have jets di 
beamed IC/CMB X-rays, which require a fast jet close to the line of sight (e.g., 
from differences in intrinsic jet power, and these effects are not easy to separate. 



broadly in accordance with recent studies (e.g., VBS03) which 
imply that the disk/corona components of RQQs do not sig- 
nificantly evolve with redshift in the range z ~ 0-5. Appar- 
ently, RLQ jets are also almost completely built-up systems at 
z»4-5. 

3.4. FIRST 0747+2739 

Although the z = 4.11 quasar FIRST 0747+2739 is radio- 
quiet (a faint detection in the FIRST survey gives a radio- 
loudness of R « 2), it is one of the most optically luminous 
quasars at z > 4 with Mb = —29.2 (compare with Figure 1 
of Vignali et al. 2003a). Richards et al. (2002) found an 
overabundance of C IV absorption systems in the spectrum of 
FIRST 0747+2739, with at least 14 independent absorption 
systems spanning a redshift range of Az « 1 longward of the 
Lya forest. They concluded that the absorption is most likely 
intrinsic to the quasar, resulting from multiple high-velocity ab- 
sorption systems at the spatial, density, or temporal edge of a 
BAL region. FIRST 0747+2739 may be a quasar in the transi- 
tion region between BAL and normal quasars. 

From the Chandra observation of FIRST 0747+2739, 
we measure an observed-frame 0.5-2 keV flux of 
1.02 x 10~ 14 erg cm~ 2 s _1 , from which we derive a rest-frame 
2-10 keV luminosity of 2.3 x 10 45 ergs -1 and an optical-to- 
X-ray spectral index of a ox = — 1 .86. Based on a significant cor- 
relation between C IV absorption and X-ray-weakness, Brandt, 
Laor, & Wills (2000) have suggested that X-ray absorption 
associated with a BAL-type outflow is the primary cause of 
X-ray-weakness in z < 0.5 blue AGN, and recent results sug- 
gest that this is likely the case for z > 4 quasars as well (see 
Vignali et al. 2003b and references therein). Despite its signifi- 
cant C IV absorption, FIRST 0747+2739 does not appear to be 
notably X-ray weak when compared to other highly luminous 
quasars (see, e.g., Figure 3 of Brandt, Schneider, & Vignali 
2004). However, since we are sampling highly penetrating X- 
rays with rest-frame energies of m 3^10 keV (and have a lim- 
ited number of photons), the X-ray column density constraints 
are not tight. Spectral fitting as described in §3.2 indicates 
that intrinsic neutral absorption with Nh < 5.5 x 10 22 cm~ 2 is 
consistent with our data. Of course, even larger column den- 
sities can be accommodated for ionized X-ray absorption (as 
is likely expected from the observation of C IV absorption). 
Thus, our results are at least consistent with expectations from 
Brandt et al. (2000), even if they cannot directly support these 
expectations. 

4. SUMMARY AND FUTURE WORK 

We have presented new Chandra detections of six RLQs and 
one RQQ at z~ 4.1-4.4, thereby significantly enlarging the 
number of z > 4 RLQs with X-ray detections. Our RLQ tar- 
gets have moderate R parameters of « 40-400; they are more 
representative of the overall RLQ population than many of the 
high-/? (R > 500) RLQs at z > 4 previously studied in X-rays. 
Our main results are the following: 

1 . Moderate-/? RLQs at z > 4 have enhanced X-ray emis- 
sion relative to RQQs of the same redshift and optical 
luminosity. The degree of enhancement, a factor of ~ 2 
in X-ray luminosity, is consistent with that observed 

d along a larger angle to the line of sight, thus making it less likely to observe 
is & Krawczynski 2002). However, the lower radio loudnesses may also result 
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at low redshift for comparable objects and is presum- 
ably associated with a jet-linked emission component. 
As in the local universe, the X-ray-to-optical flux ratios 
of moderate-/? RLQs at z > 4 are intermediate between 
those of RQQs and those of high-/? RLQs. 

2. Joint X-ray spectral analyses of moderate-/? RLQs show 
that, on average, their spectra can be fitted with a 
power-law model and Galactic absorption; the aver- 
age best-fit power-law photon index of Y = 1 .65 ±0.15 
is consistent with that for comparable objects at low 
redshift. Adding an absorption component to the fit, 
we find tentative evidence (90% significance) for ex- 
cess intrinsic X-ray absorption on average, and our 
best-fit estimate of the cold absorption column of 
Nh = 2.4+j g x 10 22 cm~ 2 is consistent with previous 
claims of increased absorption toward high-redshift 
RLQs. The somewhat higher best-fit photon index of 
r = 1.801q 25 m tn ^ s m °del is still consistent with pre- 
vious estimates for RLQs at low redshift. 

3. We do not detect extended X-ray emission associated 
with jets from any of our RLQs, which, if dominated 
by the IC/CMB process as suggested, should be signifi- 
cantly enhanced at z > 4. If the jets in these objects are 
similar to the more modest X-ray jet recently discovered 
from the z = 4.3 blazar GB 1508+5714 (which accounts 
for only as 3% of the quasar X-ray luminosity), we ex- 
pect that the jets are below the sensitivity limit of our 
snapshot observations. 



4. Overall, our data are consistent with the idea that the 
dominant X-ray production mechanisms of z > 4 and 
z ~ RLQs are similar. 

5. Richards et al. (2002) suggested that an abundance of 
C IV absorption systems in the optical spectrum of 
FIRST 0747+2739 are intrinsic to the quasar and that 
it represents a class of objects in a transition phase be- 
tween normal and BAL quasars. From our observations, 
FIRST 0747+2739 does not appear to be notably X-ray 
weak, but we cannot rule out neutral intrinsic absorption 
columns with Nh < 5.5 x 10 22 cirr 2 . 

Finally, several of the RLQs we have detected with Chandra 
are among the X-ray brightest non-blazars known at z > 4 (see 
Figure 4a). These are prime targets for follow-up X-ray spec- 
troscopy with XMM-Newton (with the caution that in some 
cases nearby bright sources may contaminate the results). Such 
spectroscopy can reveal if the majority population of RLQs at 
high redshift shows intrinsic X-ray absorption. 
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State President's Fund for Research (LCB), NASA LTSA grant 
NAG5- 13035 (WNB, DPS), and NSF grant 03-007582 (DPS) is 
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the Italian Space Agency under contract ASI I/R/057/02. 
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